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Abstract 


We derive a deterministic mathematical model that scrutinizes the dy- 
namics of cholera pathogen carriers and the hygiene consciousness of in- 
dividuals, before the illness, during its prevalence, and after the disease’s 
outbreaks. The dynamics can effectively help in curtailing the disease, but 
its effects had less coverage in the literature. Boundedness of the solu- 
tion of the model, its existence, and uniqueness are ascertained. Effects 
of cholera pathogen carriers and hygiene consciousness of individuals in 
controlling the disease or allowing its further spread are analyzed. The 
d 
d 


erential transformation method is used to obtain series solutions of the 


mT 


oe 


erential equations that make the system of the model. Simulations of 


*Corresponding author 


Received 24 May 2023; revised 31 January 2024; accepted 14 February 2024 
Sani Fakai Abubakar 
Department of Mathematics, Faculty of Physical Sciences, Kebbi State University of 


Science and Technology, Aliero, Nigeria. saniabufakai@hotmail.com 


Mohammed Olanrewaju Ibrahim 
Department of Mathematics, Faculty of Physical Sciences, University of Ilorin, Horin, 


Nigeria, moibraheem@yahoo.com 


How to cite this article 

Abubakar, S.F. and Ibrahim, M.O., Dynamics of Cholera Pathogen Carriers and Effect 
of Hygiene Consciousness in Cholera Outbreaks. Iran. J. Numer. Anal. Optim., 2024; 
14(2): 449-474. https://doi.org/10.22067/ijnao.2024.82434.1269 


449 


Abubakar and Ibrahim 450 


the series solutions of the model are carried out and displayed in graphs. 
The dynamics of the concerned state variables and parameters in the model 
are interpreted via the obtained graphs. It is observed that higher hygiene 
consciousness of individuals can drastically reduce catching cholera disease 
at onset and further spread of its infections in the population, this in turn, 


shortens the period of cholera epidemic. 


AMS subject classifications (2020): Primary 92D30; Secondary 92B05. 


Keywords: Differential transform; carriers; simulation; hygiene conscious- 


ness; series. 


1 Introduction 


Mathematical modeling can be defined as the process of selecting and utilizing 
relevant mathematical ideas to scrutinize and portray situations in a better 
way to understand and apply safer decisions about physical phenomena, and 
the entire real-life problems. Models are believed to aid in understanding var- 
ious ways phenomena and real-life problems operate. This includes epidemic 
processes, and it helps in designing effective control strategies [1]. 

The world has witnessed a remarkable improvement in hygiene, sanitary 
measures, and effective treatment of infectious diseases. This has drastically 
reduced the threat of diseases like cholera in many countries. In addition 
to this, many studies, experiments, and theories that were carried out and 
developed now and in the past have aided in reducing cholera mobility [8]. 
A report from “acap” stated that “there have been 98,9382 suspected cholera 
cases, that caused 3,293 deaths from these cases in Nigeria and that 5,343 
confirmed cholera cases with 156 deaths in Niger Republic as of October 
2021”. 


Due to poor sanitation systems, the shortage of good drinking water, 
and the high density of the population, Africa become a victim of cholera 
outbreaks, leading to a higher risk of outbreaks, and the spread of the disease. 
This coupled with shifting priority to COVID-19 resulted in the diversion of 


most resources and attention to it, leaving other health problems. This in- 
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turn caused lapses in the entire health sector and affected the process of 


contact tracing and thorough management of cholera [27]. 


1.1 Cholera 


Cholera is an illness caused by gram-negative, nearly comma-shaped bacillus 
bacteria called Vibro cholerae. The sickness’s symptoms are mostly, severe 
gray color diarrhea, with flecks of mucus, extreme thirst, drop in blood pres- 
sure, abdominal pain, and kidney failure [29]. It usually comes with vomiting 
and can dehydrates the patient in quite short period of time if prompt treat- 
ment measures are not taken. Its outbreaks usually comes after natural and 
artificial catastrophes that disrupt the public strong sanitation service sys- 
tems. The disease is still a threat and problem to global public health [26], 
especially in developing countries. 

Mishandling faeces and vomit of infected individuals easily trigger spread 
of the disease. This is so because Vibro cholerae bacteria that pass through 
the human digestive system are more infective and stubborn to control by 
human body immunity than others. Among its crucial characteristics is that 
it can survive even with limited oxygen and iron [17]. It has been estimated 
that (though, depending on the individual’s age, stomach acidity and blood 
group), an individual has to ingest up to a dose range of 10° to 10° Vibro 
cholerae bacteria cells, before developing cholera disease. The disease’s in- 
cubation period ranges (that is, from ingestion of Vibro cholerae bacteria to 
cholera infection) from 12 hours to 5 days [11, 20]. The range is pegged at, 
from 18 hours to 5 days, in [24]’s study (as cited in [12]). Vibro cholerae’s 
mechanism of causing cholera was described by [25] that usually start with 
surviving the stomach acidity, shedding its flagellum, swim deep into in- 
testine mucus. It colonizes and secretes a toxin that stimulates intestine 
mucosal cells to pump body water due to osmosis. This leads to dehydrating 
the patient. 

The onset outbreak of cholera is usually caused by the consumption of 
contaminated water. The subsequent spread of the disease can be through 


contaminated food, contact with cholera patients and unhygienic handling of 
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faeces, and vomit of cholera patients. Most mathematical models on cholera 
we visited, have not categorically scrutinized cholera pathogens since the 
pathogens cannot on their own spread without human contribution. It is 
hence vital to study the dynamics of Vibro cholerae as causative agents or 
pathogens that cause cholera disease [28]. 

Cholera is an extremely virulent disease that affects all ages and economic 
classes of a nation. It has been estimated that about two billion people 
drink water from sources that may be contaminated with faeces and that 
1.3 billion people are at risk of cholera disease in cholera endemic areas and 
regions of the world [29]. 75% of cholera pathogens carriers, manifest no 
symptoms but spread the pathogen to their immediate environment through 
faeces for a period of 7-14 days, of ingesting the pathogens by the carriers 
[13]. People with low gastric acidity are at high risk of catching cholera. 
Therefore transmission of the disease in such a class of people is rampant 
[1]. The study intends to apply mathematical modeling to best portray the 
transmission thresholds of the infections by Vibro cholerae in the host, and 


the magnitudes it attains so as to get the best ways to eradicate it. 


1.2 Related literature 


Elhia et al. [10] proposed a mathematical model on the global stability of 
the SIR cholera epidemic model that dealt with two infectious stages and 
treatment of cholera. They carried out an analysis of the model’s global 
stability, obtained Ro, and illustrated the analytical results numerically. 

A mathematical model on the controls of cholera, including hygiene con- 
sciousness as a control strategy and a compartment of the model was formu- 
lated in [21]. 

Mathematical models on the dynamics of individuals who ingested cholera 
pathogens but have not developed any symptoms of cholera disease are very 
rare. These categories of individuals unknowingly spread the pathogens of 
the disease in the environment. This in turn causes the spread of the dis- 
ease. Therefore, incorporating all possible and necessary compartments in a 


mathematical model, will surely warrant a chance to determine the essential 
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dynamics it entails. This will also help in determining the basic mechanisms 
of transmission of the disease. Its outbreak and ways individuals can be in- 
fected by the disease. It can also help in exploring various ways of its spread, 


methods of controlling it, and what can make it escalate or persist. 


These motivate us to formulate a deterministic mathematical model of 
cholera by marrying [21]’s SICRB cholera model with [10]’s SIUTR cholera 
model, to come up with SICUTRB deterministic mathematical model of 
cholera, in order to investigate the dynamics, spread, and control measures of 
cholera disease, especially through hygiene consciousness. The model has two 
incidences, linear and nonlinear. The linear incidence represents the human- 
to-human transmission of cholera, and the nonlinear is the logistic response 


_B_ 


to the increase in the cholera pathogen B, and ¢ is represented by vl gyq. 


Differential transform method (DTM) is a numerical-analytical (others 
call it semi-analytical) technique. It is used in solving linear and nonlinear, 
ordinary and partial differential equations, usually with known initial values. 
It is carried out by obtaining the approximate series solutions of the equa- 
tions. It was first used by [30] in solving linear and nonlinear initial value 
problems in electrical circuit analysis. Pukhov earlier proposed differential 
transform by computing the image of a transformed function using differen- 
tial operations, in his work ’Differential Transformations of Functions and 
Equations’ published in 1980. Zhou pioneered using the method to solve the 


initial value problem [5]. 


The method is an iterative procedure for getting analytic Taylor series 
solutions of differential equations by iterations [4, 31]. In other words, DTM 
solves an induced recursive equation from a given differential equation by 
determining coefficients of the Taylor series of the concern function [2]. Its 
analytical solutions can be constructed in the form of polynomials, whose 
obtained solutions can be presented in terms of convergence series. This in 
turn makes it easier to compute especially nonlinear equations. Application of 
the method reduces the bulk computation and helps arrive at more accurate 
approximate answers while solving linear and nonlinear differential equations. 
In recent years, the method has been employed in solving one-dimensional 


planar Bratu’s problem, differential-difference equations, delay differential 
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equations, differential-algebraic equations, and integro-differential systems 
[32]. 


Some basic definitions and properties of DTM as obtained in related lit- 
erature, such as [2, 3, 14, 18, 19, 23], are as follows: The transformation of 
the bth derivative of a function with one variable is 


b 
Gib] = ila (1) 


Here g(t) is the original function, while G[b] is the transformed function, 


and the inverse transformation is given as 


g(t) = 5 G[b|(t — to)”. (2) 


Some functions and their respective differential transformations are dis- 


played in Table 1. 


2 Model formulation 


State variables and parameters used in the model are described in Tables 2 


and 3, respectively. 
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Table 1: Some functions and their respective differential transforms 


Original Function Transformed Function 
gt)=f()+h(t) Gh =F) +A] 


g(t) = mf (#) G[b] = mF |b] m is a constant 
g(t) =1 G[b] = 6[6] 
g(t) =t G[b] = 6[b — 1 
b 
g(t) = f(t)h(t) Gb] = 5° FlalH[b - a] 
iv 
g(t) =e** OO ar 
gi =e" Gll=5 
ja) G[b] = (6+ 1)G[b+ IY 
raiO 
g(t) = °7 G[b] = (b+ 1)(b+ 2)G[b + 2] 
a 1, ifb=a 
g(t) =t Gt = 40-0) = {Te 


where 6 is the Kronecker delta 


g(t) = (1+2)° oy = cD eb) 


b! 
7 
g(t) =sin(Ot+a) G[b] = 7 sin (3 + a) 
b 7 
g(t) = cos(@t +a) G[b] = 77008 (5 + a) 


Table 2: State variables description 


Variable Description 


Susceptible Individuals 

Infected Individuals 

Hygiene Conscious Individuals 
Untreated Vibro cholerae Carries 
Infected Individuals Under Treatment 


Recovered Individuals 


wmayeiaqagan 


Concentration of Vibro cholerae Bacteria in Food 


and Water Consume by the Human Population 
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State Variables Description Value Reference 
& Parameter 
S Susceptible Individuals 900 Estimated 
Infected Individuals 50 Calculated 
Cc Hygiene Conscious Individuals 500 Estimated 
U Untreated Vibro cholerae Carries 25 Calculated 
db Infected Individuals Under Treatment 22 Estimatec 
R Recovered Individuals 3 Calculated 
B Concentration of Vibro choleraeBacteria in Food 600 Estimatec 
q Recruitment rate of human population 0.104 Estimatec 
6 Per capital natural human death rate 0.0101 [10 
U1 Rate of ingesting Vibro Cholerae bacteria from 0.5 21 
contaminated food and water consume 
by the human population 
U2 Rate of getting Vibro cholerae bacteria 03 21 
through human to human interaction 
0 Recovery rate of hygiene individuals 0.38 21 
h Hygiene consciousness of individuals 6 Estimated 
B Rate of Vibro cholerae decrease .78 Estimated 
n Rate at which infected individuals increase 10 6] 
Vibro cholerae bacteria in the environment 
f Rate of freedom from Vibro cholerae 7 Estimated 
by the hygiene conscious individuals 
k Concentration of cholera pathogen that yields 50% 10° [6] 
chance of individual developing cholera disease 
n Progression of recovery of treated individuals 0.08 Estimated 
to untreated Vibro cholerae carriers 
g2 Rate of interruption of treatment 0.2 [10] 
ry Recovery rate of untreated Vibro Cholerae carriers 0.8 Estimatec 
r2 Recovery of treated cholera patients 0.9 Estimatec 
€ Disease induced death rate for individuals in J class 0.105 —Estimatec 
LL Disease induced death rate for individuals in U class 0.1 Estimatec 
a Disease induced death rate for individuals in T class 0.08 Estimatec 
m Rate of movement of the infected class 0.4 [10] 
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(Lam \. 


Figure 1: Schematic diagram of the model 


2.1 Basic assumptions of the model 


The following basic assumptions are used in formulating the model. They 


serve as conditions under which the model will be meaningful and work effi- 


ciently. 


. A homogeneous human population is used. 


. Infected individuals are classified into under-treatment and untreated 


healthy Vibro cholerae carriers. 


. State variables and parameters of the model are positive and remain 


positive. 


. Only few cholera pathogen carriers develop severe cholera disease 


(about 75% of cholera infected individuals are healthy cholera pathogen 


carriers) during a cholera outbreak. 


. Recovery of severe cholera patients is through treatment and asymp- 


tomatic cholera pathogen carriers can heal untreated. 


. Recruitment into the population is through birth and immigration while 


the exit is through migration and natural and disease-related deaths. 


These assumptions were used to draw the flow diagram in Figure 1, then 


used to derive the model’s system of differential equations, as follows: 
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ds 


4 u1oS — (1— h)veST — 6S, (5) 
dl 
TT UGS + (1 — AjoeST —(S+etritra+ hi, (6) 
dC 
Bu e+ F +e, (7) 
dU 
pe Aa Ti+ OF AC + get — (Str tre t+ gu, (8) 
dT 
BTM tn —O+at go+ra)t, (9) 
dR 
pe it ral + rel + (ri + 172)U — oR, (10) 
dB 
di = (1 h)nl BB. (11) 


2.2 Existence and uniqueness of solution of the model 


Formulated mathematical models need to undergo a surety test to ascertain 
the existence and uniqueness of their solution. Using Lipschitz criteria as in 
the definition of basic terms, the existence and uniqueness of a solution of 
the model are established as follows: 

Theorem 1 stated below, is from Derrick and Grossman (1981), also available 
in [9, 16], together with Theorem 2 are adopted to prove the existence and 


uniqueness of solution of the model. 


Theorem 1. Suppose that D! denotes the region 


|t—to |<, || a — ao ||< m, 
@ = (@1,42,...,@n), ao = (410, @20,---,4no); 
and that 
Ilg(a@1,t) — g(aa,t)|| < plla1 — al). (12) 


If the pair (a;,t) and (az,t) belong to D!, where p is a positive constant, 
then there is a constant 6 > O such that there exist a unique continuous 
vector solution a(t) of the system in the interval t — to < 6. This condition 


is satisfied provided that 
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Og: ; 
E. (a bee eee 
] 


is continuous and bounded in D!, (where i represents equations in the system 


and j represents the compartments). 


Now, considering the region 0 < w < R. Bounded solution of the system 
of equations (5)—(11) will be looked for in the same region whose partial 
derivatives satisfy 6 < w < R, where 6 and w are positive constants, as in 
the theorem that follows. 


Theorem 2. Let D denotes the region 0 < w < R. Then the system of 
equations (5)—(11) has a unique solution if it is established that 


al i= 1,2,3,4,5, 6, 7, 
Ot; 
are continuous and bounded in D. 


The proof of the existence and uniqueness of solution of model two system 


of equations (5)—(11) are as follows: 
Proof. Let 
ly = 7 — v10~S' == (1 = h)voST = 6S, 


lg = 0, 6S + (1—h)veST—(d+e4+r +192 +h)I, 
ee eee we ce) 


y= (L—r)mI + (64+ f)C+ mT —(64+ritratmtayu, 
ls =7mI 4+ qmU —(d+a4+ 924+ 1a)T, 

Ie = (ri tre) + reT + (ri + 1r2)U — dR, 

ly =(1—h)nI — BB. 


Partial derivatives with respect to each state variables of the system (5)— 


(11) are obtained as follows: 


For ly, 
Ol, | Ol, | 

Ee |—u1@ — (1 — h)vel — 6| < on, orl = |—(1 — h)v2S] < oo, 
Oly Oly Ol, Oly Oly oe 
OC OU OT OR OB 
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Similarly, for lg, 13, l4, ls, lg, and l7, all have the absolute values of their 


partial derivatives less than infinity. 


Solution of the system of equations (5)-(11) exists and is unique in R’ 


since all the partial derivatives are bounded and defined. 


2.3 Boundedness of solutions of the model 


Let Py be the human population of the model represented by equations 


(5)-(10). The variables and parameters used are nonnegative. 


Theorem 3. All solutions of equations (5)—(10) of the model are bounded 


in R&, and solution to equation (11) is bounded in R,. 


Proof. Let 
Pe=—S+I+C+U4T +h. (13) 


Differentiating (13) with respect to t gives 


dPy dS di dC dU dv dR 
i a a ae ae ee a) 


Simplifying (14) gives 


dPy 


ape 6S —(e—-m+6)I — (6)C —(u+6)U —-(a+6)T-6R. (15) 


Taking and using any p > 0, multiplying the p by (13), then adding the result 


to respective sides of (15), and simplifying, thus we have 


o(Px) = p(S+I1+C+U+T+R), 


dP 

+ oP =9 G=ps—eHamo—p l= —peCa(to=pe 
—(a+6—-p)T — (6—p)R. 

Let 


Aq = (6—p)S—(e—m+4d—p)I—(6—p)C—(u+6—p)U—(a+d—p)T—(5—p)R 


Iran. J. Numer. Anal. Optim., Vol. 14, No. 2, 2024, pp 449-474 


Dynamics of Cholera Pathogen Carriers and Effect of Hygiene ... 


with 
p=min{d,e—m+6,6,u+6,a+ 4,6}, 
which means that p is the minimum value of 6. Hence 6 — p < 0. Similarly, 


p is the minimum values of ¢ —-m-+ 6, 4 +6, and a+ 06. This implies 
e—-m+6—p<0,u+d-—p<0,anda+d—-—p<0. 


Therefore ; 
Pr 
—+pPy=q-A 
rr: + pF = q q 
but A, <0. Taking Ag = 0 gives 
dP. 
at PH S 4. (16) 


Now, using the differential inequality theorem from Gronwall (1919) and 
the Brikhoff and Rota version as used in [7, 15, 21, 22], (16) gives 


eee 17 
He 20 Oke (17) 


Separating the differential inequality (17) gives 


dPy 
—=— < dt. 18 
q —pPH ~ oe 
Integrating the both sides of (18) gives 
1 
ae (q—pPyH) <tte. (19) 
Solving (19) in terms of c gives 
1 
Sing ge Se (20) 
p 


Evaluating (20) at = 0 and Py = Py(0) gives 
1 
“a (q — pPx(0)) <e. (21) 
Substituting (21) into (19) gives 
1 1 
—2in(q— pPu) <t+ (—Lim(q—pPa(0))) (22) 


Multiplying (22) by —p gives 
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In (q— pPx) < —pt + n(q— pPx(0)). 


Taking exponent of (23) gives 


q—pPu < e” (q— pPx(0)) 


ast ow, 


ef! +0. 


Also, 


e* (q— pPx(0)) + 0, 


q—pPy <0, 


q 


462 


(23) 


(24) 


so Py < me All state variables and parameters of the model are positive. 


Hence the sum Py is greater than zero. Therefore 


pL Pes: 


This cumulatively implies the system of equations (5)—(10) is bounded and 


its solutions enter into the region 


Ah = {(S:1C.0,T.R) ERE 0< Pe < At 


(25) 


Hence, the solution of equations (5)—(10) is bounded and belongs to feasible 


region H,. On the other hand, the nonhuman population is described by 


equation (11). Its boundedness is sought as follows. 


Let Pg denote the number of cholera pathogens, so that equation (11) 


can be written as 


dt 


= (1—h)nI — BPp. 


From (13) wn <q —pPy. Since Py < then from the same (13), it 
q 
A 


implies that I < 


Therefore 


holds and can be written as 


dPz 


< dt. 
obo) — BPs 


(26) 
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Integrating the both sides of (26) gives 


sin (fa h) Pe) Kt4- ci, (27) 
Evaluating (27) at t = 0 and Pg = Pp(0) gives 
5m (“a h) 5Pa(0)) 2 (28) 


Substituting (28) into (27) gives 


5n(@a h) BP») <1 jn(@a h) 3P9(0)). (29) 


Multiplying both sides of (29) by —6 gives 


In (fa iim Pe) <—ft+In (Za hy = Pa(0)) . (30) 
Taking exponent of both sides of (30) gives 
Tl — h) — BPs < e Pt FAH) BP (0), 
p 
which can be written as 


1h) — BPp < e# (“a 3i= aP9(0)) 
p p 


Then, as t —> 00, e~F* —+ 0, and also, e~** (2a —h)- 6Pa(0)) => 0. 
it is obtained that 


qn 
Pax SH). 
BS Gol ) 


SP hs, 
Bp 
Therefore 
qn) 
0< Pp < —(1—-h). 
Bp\ 
The nonhuman population equation solution is hence bounded and enters the 
region 
= {BER,.0< Po < 2am}. (31) 
Bp 
Altogether, the feasible solutions of the system of equations (5)-(11) are 


bounded and enter the region 
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A= {(S.GU.7.R) eR{,BeER,,S,1,C,U,T,R,B> 0, 


qd q") 
0< Py <-,0<P, < fan}. 32 
eS B Bp | ) (32) 


3 Results of the model using DTM 


Gi =d+etrnt+reth, 
Gg=dO+tpt+rnt+ret+ 1, 
G3 =d+at+g.+?Po, 
Gy=o+0+f, 
G;=0+f, 

Ge =1— nh, 
G7 = 19, 

Gg = v2 — voh, 
Go = m7, 

Gio =Mm— mr. 


The system of ordinary differential equations of the model is now written as 


att = q — G7s(t) — Ggs(t)i(t) — ds(t), (33) 
a = G7s(t) + Ggs(t)i(t) — Gii(t), 84) 
att) = hi(t) — Gac(t), a) 

a = Gioi(t) + Gse(t) + gati, (t) — Gout) (38) 
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-~ eRe MT CCN Or (37) 
a) = Gsc(t) + rpu(t) + reti(t) — dr(t), (38) 
ee ee 2?) 


dt 


Let S(b), I(b), C(b), U(b), T(b), R(b), and B(b) denote the differential 
transformation of s(t), a(t), c(t), u(t), ti(t), r(t), and 61(t), respectively. Ap- 
plying the fundamental operations of DTM provided in Table 1, the following 


recurrence relations are obtained: 


3 1 a=0 
b 
I(b+1)= a (ase + Gg S° (S(a)I(b- a)) - curt) , (41) 
a=0 
1 

C+ 1) = FH (hI) - GsC(), (42) 
U()+1) =; (Gol) + GsC() + 927) - GV), (48) 
T(b+1) = 55 (Gollb) + nV) - GsTO)), (44) 
ROV+1) = p> (GsCO) +nUW) + rT) —45RO)), (45) 
Bb+1) = = Goll) - 6B). (46) 


Substituting the initial conditions 

s(0) = 900, i(0) = 50, c(0) = 500, u(0) = 25, t1(0) = 22, r(0) = 3, b1(0) = 
600 

and the values of the parameters 

q = 0.1, 6 = 0.0101, vy = 0.5, vg = 0.03, 0 = 0.38, h = 0.6, 8 = 0.8, 7 = 10, 
f = 0.7, g: = 0.08, g2 = 0.2, 71 = 0.8, ro = 0.9, € = 0.1, w = 0.1, a = 0.08, 
m = 0.4, rT = 0.5, d = 0.000599, 

into DTM conditions (40) to (46). Using Maple 18, we obtain 
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S(1) = —549.25955, ——«‘(1) = 419.76455, 

C(1) = —515.05, U(1) = 507.1475, 

$(2) = —2099.094692, (2) = 1596.031182, 

C(2) = 406.6573675,  U(2) = —716.84851, 

5(3) = —4396.37805,  1(3) = 3121.246753, 

C(3) = 171.4405043, -U(3) = 709.150619, 

S(4) = —2483.196769, (4) = 613.668424, 

C(4) = 421.4651895,  U(4) = —129.7781802, 

S(5) = 15561.71779, (5) = —15852.50219, 

C(5) = —18.24762974, U(5) = 170.74677881, 

S(6) = 56847.30013, (6) = —50505.80978, 

C(6) = —1581.934929, U(6) = —585.9512738, 

T(1) = -14.1822, R(1) = 579.7697, B(1) = 320, 

T(2) =70.7014731, -R(2) = —84.577827,  B(2) = 3230.1164, 

T(3) = 59.23884417, _R(3) = —23.26776309, B(3) = 7650.80193, 

T(4) = 152.6203129, R(4) = 201.506551,  —-B(4) = 1095.82663, 

T(5) = —13.8564008, R(5) = 97.33658518,  _B(5) = 210.966696, 

T(6) = —523.3916988, R(6) = 17.2393537, _B(6) = —42301.46807, 
R(7) = —378.3539126, B(7) = —110607.3974, 


The series solutions of the transformed equations b ranges from 1 to 6 or 7: 


s(t) = 5_ t°S(b) = 900 — 549.25955¢ — 2099.094692¢? — 4396.37805¢° 
— 2483.196769t* + 15561.71779¢° + 56847.30013¢° + --- , 
i(t) = S_ t'1(b) = 50 + 419.76455t + 1596.031182¢? + 3121.246753¢" 
+ 613.668424t* — 15852.50219t° — 50505.80978¢° + --- , 
j= S- t°C(b) = 500 — 515.05t + 406.6573675t? + 171.4405043¢3 
+ 421.4651895t* — 18.24762974t° — 1581.934929t° +.--- , 
i)= $ t?U(b) = 25 + 507.1475t — 716.84851t? + 709.150619¢3 


— 129.7781802t* + 170.74677881t° — 585.9512738¢° + 
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b 
ti(t) = > t?T(b) = 22 — 14.1822t + 70.7014731t? + 59.23884417¢° 
a=0 
+ 152.6203129¢* — 13.8564008¢° — 523.3916988¢° + --- , 
b 
r(t) = > t?R(b) = 3 + 579.7697t — 84.577827t7 — 23.26776309¢° 
a=0 


+ 201.506551t* + 97.33658518¢° + 17.2393537¢° 
— 378.3539126t7 + ---, 


7 
bi(t) =S_ tB(b) = 600 + 330¢ + 3226.6364t? + 7662.06045¢° 
b=0 


+ 10960.11186¢* + 214.6803636t° — 42259.50383t° 
= 110607.3974¢' +++. 


The Maple 18 plots of s(€), i(¢), c(t), u(t), t1(t), r(t), and b,(t) are displayed 
in Figures 2-6. 


4 Conclusion 


The model proposed how hygiene consciousness of individuals influences the 
susceptible and infected individuals. A high rate of h, which is hygiene con- 
sciousness of individuals, increases number of susceptible individuals while 
the low rate of h decreases their number, for details, see Figures 3, 4 and 
5. The untreated Vibro Cholerae carriers, are seen physically as healthy, be- 
cause they do not manifest any symptoms of cholera, but spread the pathogen 
without knowing. The rate of freedom of individuals from Vibro Cholerae was 
determined by their hygiene consciousness (see Figure 6). In Figure 3, the 
number of susceptible individuals reduced from 900 to about 100 individuals 
in less than a week, infection reduced to zero in about half of a week, the 
peak of the infection is 350 out of 900 susceptible individuals when the rate of 
hygiene consciousness is 63%. While the number of susceptible individuals re- 
duced from 900 to 90 individuals, in about two weeks’ time, infection reduced 
to zero in two weeks, while the peak of infection is 200 out of 900 susceptible 


individuals, when the rate of hygiene consciousness is 82%. The number of 
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susceptible individuals dropped from 900 to 400 individuals, in about a week, 
infection reduced from 50 to zero individuals in the same period, peak of the 
infection is 50 out of 900 susceptible, when the rate of hygiene consciousness 
is 99%, as in Figure 5. Reduction of hygiene consciousness results in the pro- 
longation of the cholera illness period to about two weeks and the maximum 


number of infected individuals reach 200. 


1600 


3 8 


S(H.10.CO.UO.TOR & BO 
EH 


0 


0 0s 1 15 25 
Time in Weeks 


— Susceptible — Infected 
Hygiene Conscious Untreated V. cholerae Cariers 

— Infected Under Treatment —— Recovered 

| V. cholerae Population 


Figure 2: Combined dynamics of the seven compartments of the model plotted against 


time. 


Dynamics of Hygiene Consciousness, at h = 0.63 


000 


Susceptible Infected & Hygiene Conscious Individuals 


0 01 02 06 07 08 


03 04 OS 
Time in Weeks 


Susceptible —— Infected —— Hygiene Conscious 


Figure 3: Dynamics of hygiene consciousness, its effect on susceptible and infected 
individuals at h = 0.63. As higher hygiene consciousness is maintained, number of 


individuals susceptible to cholera drastically reduced, and cholera infection stopped. 
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Dynamics of Hygiene Consciousness, at h = 0.82 


Susceptible Infected & Hygiene Conscious Individuals 


0 05 1 1S, 2 
Time in Weeks 


Susceptibles ——— Infected ——— Hygiene Conscious 


Figure 4: Dynamics of hygiene consciousness, its effect on susceptible and infected 
individuals at h = 0.82. Maintaining higher rate of hygiene consciousness, will make the 
number of individual susceptible to cholera reduced drastically and cholera spread and 


infection stopped immediately. 
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Dynamics of Hygiene Consciousness at h=0.99 
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Figure 5: Dynamics of hygiene consciousness, its effect on susceptible and infected 
individuals at h = 0.99. Applying and maintaining about 100% hygiene consciousness, 
in societies, will within just about 10 days, make the number of susceptible individuals 


to approach zero as well as the rate of cholera infection. 
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Figure 6: Effect of hygiene consciousness on untreated Vibrio cholerae bacteria carriers. 
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